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We investigate the non-Abelian Bianchi identities in pure SU(2) lattice Yang-Mills theory in three 
and four dimensions. The non-Abelian Stokes theorem proposed recently allows to formulate the 
Bianchi identities in terms of local physical fluxes. Then the violation of Bianchi identities be- 
comes a well defined concept ultimately related to chromomagnetic fields degeneracy points. 
We present numerical evidences that in D=4 the suppression of the Bianchi identities violation 
destroys confinement while the removal of the degeneracy points drives the theory to the topolog- 
ically non-trivial sector. 
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Bianchi Identities and Chromomagnetic Fields Degeneracy Points on the Lattice 

The possibility of Bianchi identities violation had always been considered as probable cause of 
confinement (see, e.g., Ref. ^ for recent discussion). However, detailed treatment of the problem 
in non-Abelian case was missing until recently [^] and here we briefly review the main points of 
this paper. Formally, the Bianchi identities in three and four dimensions 

DpPpv = 0, F« v = X -e^ Xp F a Xp (AD) , D.-S, = 0, Bf = l/2e ijk F? k (3D) (1) 

look similar, but in fact there is a crucial difference: for given non-Abelian field strength Eq. (|l|) 
allows to express the gauge potentials as local single- valued functions of F^ v , while in three di- 
mensions it does not constraint A^ at all. The inversion A(F) is possible away from the points 
chromomagnetic fields degeneracy !§], where detT = 0, Tjfi = e abc F^ v . Note that 

det T - det , £ MV = K vtl = I e ahc F« p F h pX F c lv (2) 

and each element of K^ v is gauge invariant determinant constructed from chromoelectric Ef and 
chromomagnetic B" fields. It is amusing that our investigation of the Bianchi identities ([[]) and 
possibility of their violation naturally leads to the same determinants (||). Moreover, the degeneracy 
points turn out to be closely related to the gauge fields topology. 

We start from the observation [Q] that the Bianchi identities ([]]) are equivalent to the require- 
ment that the non-Abelian Stokes theorem (NAST) being applied to closed infinitesimal surface So 
always gives unity on the r.h.s. (surface independence) 

P s exp^ J ^ v d 2 a^ v = l=e idn - 2 ^, n 2 = 1 , q G Z, (3) 

S ,8S =0 

where for definiteness we used the operator formulation ^ of NAST and r.h.s. was represented 
in most general form. The color direction n is gauge variant and is of no concern. However, the 
possibility to give unambiguous gauge invariant meaning to the integer q ("magnetic charge") is 
important since if it is non vanishing then the continuum Bianchi identities are violated somewhere 
inside So- Note that the term "magnetic charge" here has nothing to do with its usual meaning, in 
particular, its conservation is neither assumed nor implied (see also below). Since Eq. @ is the 
definition of the magnetic charge the only way to make sense of q is to express the non-Abelian 
Stokes theorem in terms of gauge invariant physical fluxes. 

The clue to the above problem is to choose the appropriate parametrization [(J for the Wilson 
loops 

fT+t 

W(C,t)=Vexpij A(r)dT = exp{idn(C,t)-^(C)}, (4) 

n 2 (C,t) = l, W(C)= ^TrW{C,t) =cos<£(C), 

where C = {x(t),0 <t< T,x(0) = x(T)} is some closed contour and it is understood that W(C) ^ 
±1 genetically. The flux magnitude <&(C) € [0;n] is gauge invariant, while the instantaneous flux 
color direction n(C,t) explicitly depends on t and is in the adjoint representation. Note that the 
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flux magnitude is insensitive to the change of contour orientation while the flux direction changes 
sign. Consider now another contour C' which touches (or intersects) C at point x(t) = x'(t r ). Ev- 
idently, the relative orientation (angle in between) of n(C,t) and n(C',t') remains intact under the 
gauge transformations. Moreover, the construction could be iterated: for N contours C,, i = 1, ...,N 
intersecting at one point the relative orientation of instantaneous fluxes «W at that point is gauge 
invariant. In particular, the oriented area 1 ^(n^, ...,n^) of spherical N-vertex polygon is well 
defined. 

The next point concerns the relation between the physical field strength and the corresponding 
infinitesimal Wilson loop. The basic and well known fact is that the lattice area element dx^dx v 
is unoriented dx^dx y = dx v dx^ contrary to the usual continuum relation 8o^ v = dx^ f\dx v = 
— dx v A dx^- . Therefore in order to define the lattice counterpart of the field strength a canonical 
orientation of all elementary plaquettes should be fixed first. Moreover, the canonical ordering is 
well known, the conventional agreement is to consider 8o^ v with /I < v only. From now on we 
always assume that the infinitesimal fluxes are constructed with canonically oriented plaquettes. 

Once the canonical orientation is fixed, the NAST presented and discussed in details in Ref. [|]] 
becomes completely unambiguous. It relates the flux piercing the large contour C with magnitudes 
and relative orientations of elementary fluxes on (arbitrary) surface bounded by C. The requirement 
of surface independence (^) considered for elementary lattice 3-cube c becomes then 

£ £ Ha^Vf ,4 3) ) =2nq(c), (5) 

pddc x€8c 

where <&(/?) is the plaquette flux, 1 /2TrW(p) = cos Q.i, was discussed above and accounts for 
the difference in color orientations of three fluxes emanating from 3-cube and having one common 
point x £ 8c. The factors I(p) = ±1 are the usual incidence numbers: I(p) = 1 if vertices of p are 
followed in the canonical order and I(p) = — 1 otherwise. Eq. (||) is the definition of the magnetic 
charge associated with 3-cube and indeed expresses the essence of non-Abelian Bianchi identities. 
In particular, the continuum Bianchi identities ([]]) are violated inside c if q(c) 7^ 0. 

It is important that the non-Abelian Bianchi identities ([]]) became completely Abelian-like @ 
in our approach. Moreover, it is possible to construct a specific but unique cell complex in which 
every term on l.h.s. of Eq. (Q) is associated with particular 2-dimensional cell. Then Eq. (|5|) is the 
usual coboundary operation d : C 2 — > C 3 acting on particular 2-cochains. One could say that the 
non-Abelian nature of ([]]) had been traded for much more complicated geometry underlying (§), 
but which allows purely formal investigation. In particular, it is clear that the 2- and 3-skeletons of 
the cell complex are not exhausted by the original plaquettes and 3-cubes (see Ref. [§] for details). 
Then the study of d : C 2 — > C in its generality leads to the consideration of the "magnetic charges" 
assigned to all the 3-cells of the complex. It is true that some of these "new" 3-cells are trivial and 
the corresponding magnetic charge is identically zero. However, there are a non-trivial 3-cells as 
well which are distinct from the original lattice 3-cubes. One could show that the corresponding 
charges q signal the zeros of the determinants entering (Q). Formally the charges q are coming from 
the Q. terms in (||) and indeed depend only upon flux color directions, not their magnitudes. Note 
that the union of all magnetic charges q, q is indeed conserved, but the conservation takes place not 
in the original hypercubical geometry. 

1 The total area of unit two-dimensional sphere is taken to be In. 
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Figure 1: The correlation function of Polyakov lines at y = and various y in three (left) and four (right) 
dimensions. 

To summarize, the attempt to consider gauge invariant content of the Bianchi identities ([]]) 
forces us to introduce two types of the "magnetic charges" q, q. First one, Eq. (||), is responsible 
for the Bianchi identities violation and is assigned to lattice 3-cubes. Second one, being assigned to 
lattice sites, is directly related to the degeneracy points, where a particular determinants constructed 
from Ef, B" vanish. 

Numerical Experiments 

Simple weak coupling analyzes show that in perturbation theory the q, q charges form tightly 
bounded dipole pairs, but cannot annihilate because at any finite lattice spacing they by definition 
are finitely separated. This invalidates the consideration of the corresponding densities since they 
both would diverge in the continuum limit. On the other hand, the definition of q and q charges 
is local and gauge invariant. Therefore, nothing prevents us from modifying the Wilson action to 
include additional terms which suppresses the charges q, q 

S = -/3£ lTrV p + YZ\9(c)\ + ?Z\9(.')\ «> 

p c s 

and allows to study their physical relevance. The particular limit y — ► oo is of special interest since 
it corresponds to the theory with nowhere violated Bianchi identities. As far as the y coupling 
is concerned the limit y — > oo might not correspond to sensible theory. For instance, the nowhere 
vanishing detB implies that it is of the same sign everywhere, which contradicts the perturbative 
expectations and probably violates CP symmetry. 

The model (§) was studied in D = 3, 4 along the lines 7 = 0, y = at fixed value of the 
gauge coupling j3. The simulations were performed on 20 3 and 12 4 lattices at /3 = 6.0 and j8 = 2.4 
correspondingly. We implemented the Metropolis algorithm which is the only one available at 
non-zero y, y. The procedure turns out to be very time consuming. Because of this the volumes 
were not so large and only the limited set of (y, y) points were considered. At each (y, y) val- 
ues we generated about one hundred statistically independent gauge samples separated by ~ 10 
Monte Carlo sweeps. The observable of primary importance is the correlator of the Polyakov lines 
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Figure 2: Left: topological susceptibility in D=4 at y = and various y. Right: Monte Carlo history of 
global topological charge at y = and two y values. 



(P(0)P(R)) for which the standard spatial smearing [[7|] and hypercubic blocking [j|] for temporal 
links were used. In D=4 we also monitored the topological charge Q and the topological suscep- 
tibility x = (Q 2 )/V measured with overlap Dirac operator ^ (see, e.g. Ref JToj ] for review and 
further references). 

Consider the response of the heavy quark potential on rising y coupling. In three dimensions 
(Fig. [I], left) the correlator (P(0)P(R)) shows almost no sign of y coupling dependence, in par- 
ticular, the asymptotic string tension at large y is equal to its value in the pure Yang-Mills limit. 
However, the situation changes drastically in D=4. It is apparent (Fig. [lj, right) that the correlation 
function (P(0)P(R)) tends to non-zero positive value at large separations 



lim (P(0)P(tf)) r >i = const >0, 



(7) 



when 7 coupling becomes of order few units. On the other hand, other measured observables 
do not show strong dependence on y coupling. In particular, the topological charge Q stays 
at zero in average albeit with slightly narrower distribution. As is clear from Fig. || (left) the 
topological susceptibility drops at y rj 1 by approximately 25% and then stays almost constant, 
lim X 1 ^ 4 (y) = 163(8) MeV. Therefore the dynamics of YM fields in D=3 seems to be almost 

y — >oo 

insensitive to whether or not the Bianchi identities are violated. Contrary to that, in D = 4 case the 
suppression of the Bianchi identities violation seems to destroys confinement while other measured 
characteristics of the theory remain essentially unchanged. For interpretation and more details see 
Ref. [§. 

As we noted already the suppression of degeneracy points might not be physically meaningful. 
Indeed, we found that at any y > the Polyakov lines correlator becomes an oscillating function of 
the distance signaling the reflection positivity violation. In four dimensional case the natural reason 
of this is non-zero in average global topological charge. The right panel on Fig. ^ shows the Monte 
Carlo history of the topological charge on 8 4 lattice at jS = 2.30, y = 0, y = 0. 1 , 0.5 when the starting 
configuration was thermalized at y = y = (the lattice geometry and the j3 coupling were changed 
for reasons to be explained shortly). In particular, for < y 1 the mean topological charge is 
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shifted only slightly being of order few units. However, once the y coupling becomes comparable 
with unity Q flows away from zero during Monte Carlo updating towards an extremely large and 
positive values with almost constant and very high rate. In fact, it quickly becomes too large to be 
technically accessible for us and this was essentially the reason to consider so small lattices here. 
The volume dependence of (Q) could be inferred by noting that the last term in responsible for 
the rapid increase of the topological charge is the bulk quantity. Therefore (Q) is expected to be 
proportional to the volume at fixed y although we had not investigated this numerically. 

To conclude we note that the numerical results confirming the connection between confinement 
and violation of Bianchi identities should be improved in order to be convincing. In particular, we 
hope to study elsewhere the larger volumes and j3 coupling. Nevertheless, our data clearly indicates 
the relation between the degeneracy points and the gauge fields topology, which is worth to be 
further investigated. 
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